Classical versus quantum errors in quantum computation of dynamical systems 



O 

o 

(N 



> 

On 
00 

in 
o 

o 

Ok 
^— > ■ 



Davide RossinQ 

Center for Nonlinear and Complex Systems, Universita degli Studi dell'Insubria and 
NEST- INFM & Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy 

Giuliano Benent{j] and Giulio CasatO 

Center for Nonlinear and Complex Systems, Universita degli Studi 
dell'Insubria and Istituto Nazionale per la Fisica della Materia, 
Unitd di Como, Via Valleggio 11, 22100 Como, Italy 
(Dated: February 1, 2008) 

We analyze the stability of a quantum algorithm simulating the quantum dynamics of a system 
with different regimes, ranging from global chaos to integrability. We compare, in these different 
regimes, the behavior of the fidelity of quantum motion when the system's parameters are perturbed 
or when there are unitary errors in the quantum gates implementing the quantum algorithm. While 
the first kind of errors has a classical limit, the second one has no classical analogue. It is shown that, 
whereas in the first case ("classical errors") the decay of fidelity is very sensitive to the dynamical 
regime, in the second case ( "quantum errors" ) it is almost independent of the dynamical behavior of 
the simulated system. Therefore, the rich variety of behaviors found in the study of the stability of 
quantum motion under "classical" perturbations has no correspondence in the fidelity of quantum 
computation under its natural perturbations. In particular, in this latter case it is not possible 
to recover the semiclassical regime in which the fidelity decays with a rate given by the classical 
Lyapunov exponent. 

PACS numbers: 03.67.Lx, 05.45.Mt 



I. INTRODUCTION 

Fidelity is a very convenient tool to characterize the 
stability of quantum computation. It is defined as f(t) = 
\(ip(t)\ip e (t))\ 2 , where the two state vectors \ip(t)) and 
\ipe(t)) are obtained by evolving the same initial state 
\ipo), under ideal or imperfect quantum gates, respec- 
tively. Here e measures the imperfection strength and 
we assume that the perturbed gates are still unitary. If 
the fidelity is close to one, the results of the quantum 
computation are close to the ideal ones, while, if / is sig- 
nificantly smaller than one, then quantum computation 
does not provide reliable results. 

More generally the fidelity (also called Loschmidt 
echo) is a quantity of central interest in the study of 
the stability of dynamical systems under perturbations 

Biaiiiiiiiimiiiiiiiiiiaiiliiiiiiiii. 

The decay of fidelity in time exhibits a rich variety of dif- 
ferent behaviors, from Gaussian to exponential or power- 
law decay, depending, e.g., on the chaotic or integrable 
nature of the system under investigation, on the initial 
state (coherent state, mixture, etc.), and, for integrable 
systems, on the shape of the perturbation and on initial 
conditions. In particular, in the chaotic, semiclassical 
regime, and for strong enough perturbations, it has been 
shown that the decay rate is perturbation independent 
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and determined by the Kolmogorov-Sinai entropy, related 
to the Lyapunov exponent of classical chaotic dynamics 


On the other hand, the simulation of the quantum 
dynamics of models describing the evolution of com- 
plex systems promises to become the first application 
in which a quantum computer with only a few tens of 
qubits may outperform a classical computer. Indeed, ef- 
ficient quantum algorithms simulating the quantum evo- 
lution of dynamical systems like the baker's map [lflj . 
the kicked rotator [2(j, and the sawtooth map [21| have 
been found, and important physical quantities could be 
extracted from these models already with less than 10 
qubits H2, EI • Therefore, these quantum algorithms 
may constitute the ideal software for short- and medium- 
term quantum computers operating with a small number 
of qubits and the most suitable testing ground for in- 
vestigating the limits to quantum computation due to 
imperfections and decoherence effects. In this context, 
we point out that the fidelity of quantum computation 
has been evaluated for the quantum baker's map using 
a three-qubit NMR-based quantum processor j25j. We 
also note that efficient quantum algorithms to compute 
the fidelity have been proposed in Refs. [261127^ . 

From the viewpoint of computational complexity, the 
following question naturally arises: given a generic dy- 
namical system, is it possible to find its solution at time t 
efficiently, including into consideration unavoidable com- 
putational errors? We recall that the classical dynamics 
of chaotic systems is characterized by exponential sen- 
sitivity: any amount of error in determining the initial 
conditions diverges exponentially, with rate given by the 
largest Lyapunov exponent A. This means that, when 
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following a given orbit, one digit of accuracy is lost per 
suitably chosen unit of time. Therefore, to be able to 
follow one orbit up to time t accurately, we must input 
Oft) bits of information to determine initial conditions. 
On the other hand, the orbit of a non-chaotic system is 
much easier to simulate, since errors only grow linearly 
with time. Owing to the exponential instability, classi- 
cal chaotic dynamics is in practice irreversible, as shown 
by Loschmidt echo numerical simulations of Ref. [28| : 
if, starting from a given classical distribution in phase 
space, we simulate the dynamical evolution up to time 
t and then, by inverting at time t all the momenta, we 
follow the backward evolution, we do not recover the ini- 
tial distribution at time 2t. This is because any amount 
of numerical error in computer simulations rapidly ef- 
faces the memory of the initial conditions. On the con- 
trary, the same numerical simulations in the quantum 
case show that time reversibility is preserved in the pres- 
ence of small errors. 

In view of the above considerations, it is natural to 
inquire whether the degree of stability of a quantum al- 
gorithm depends on the nature (chaotic or non-chaotic) 
of the simulated dynamics. We will show that the decay 
of the fidelity of a quantum algorithm in the presence of 
perturbations in the quantum gates is almost indepen- 
dent of the dynamical behavior of the simulated system. 

In this paper, we will consider a quantum system, the 
so-called sawtooth map, which can be simulated effi- 
ciently on a quantum computer and whose underlying 
classical dynamics, depending on system's parameters, 
can be chaotic or non-chaotic. We will outline the main 
differences that occur in calculating the fidelity decay 
with "classical" and "quantum" perturbations on the dy- 
namical system. 

• By classical perturbations, we mean perturbations 
of the system's parameters that have a classical 
limit. For instance, in this paper we perturb, at 
each map step of the sawtooth model, the kicking 
strength A; by a small amount 5k(t) -C k, where 
t measures the number of map iterations. Note 
that this kind of perturbation, when applied to the 
classical motion, disturbs a given orbit by a small 
amount at each map step and therefore, to some 
extent, mimics the presence of round-off errors in a 
classical computer. 

• By quantum perturbations, we mean errors intro- 
duced at each quantum gate (in this paper, we con- 
sider unitary, memoryless errors). These quantum 
errors are unavoidable during a quantum computa- 
tion, due to the imperfect control of the quantum 
computer hardware and they do not have classi- 
cal analogue. We will show that the fidelity decay 
evaluated with quantum errors is not capable of 
distinguishing between the classically integrable or 
chaotic nature of the simulated dynamics, being es- 
sentially independent of it. 

This paper is organized as follows. In Sec.^ we briefly 



describe the sawtooth map model and a quantum algo- 
rithm which efficiently simulates it. We also introduce 
our quantum and classical error models and discuss how 
to efficiently evaluate the fidelity on a quantum com- 
puter. In Sec. 1 1 1 1 1 based on extensive numerical simu- 
lations, we analyze the differences between the fidelity 
decay in the presence of classical and quantum error. Fi- 
nally, in Sec. II VI we present our conclusions. 



II. THE PERTURBED QUANTUM SAWTOOTH 
MAP MODEL 

In order to illustrate the striking differences between 
the fidelity decays induced by classical and quantum er- 
rors, here we consider the quantum sawtooth map model. 
This map is one of the most extensively studied dynami- 
cal systems, since it exhibits a rich variety of different dy- 
namical regimes, ranging from integrability to chaos, and 
interesting physical phenomena like normal and anoma- 
lous diffusion, dynamical localization, and cantori local- 
ization eh, ill 

The sawtooth map is a periodically driven dynamical 
system, described by the Hamiltonian 

WMHy-^y^ £ ^-OT), (1) 

J= — oo 

where (n, 9) are the conjugated action-angle variables 
(0 < 9 < 27t). The time evolution t -> t + T of this 
system is classically described by the map 

n = n + k{6-ir), 9 = 9 + Tn, (2) 

where the bars denote the variables after one map iter- 
ation. By rescaling n — > p = Tn, the classical dynamics 
is seen to depend only on the parameter K = kT . The 
classical motion is stable for —4 < K < and com- 
pletely chaotic for K < —4 and K > 0: the maximum 
Lyapunov exponent is A = ln[(2 + K + V K 2 + 4K)/2] 
forX>0, A = In |(2 + K - y/K 2 + 4K)/2\ for K < -4, 
and A = for —4 < K < 0. As shown in Fig. ^ m 
the stable, quasi-integrable regime, the phase space has 
a complex structure of elliptic islands down to smaller 
and smaller scales. Note that the integrable islands are 
surrounded by a non-integrable region, and that each tra- 
jectory diffuses (anomalously) in this region. The cases 
K = 0, —1, —2, —3, —4 are integrable. 

The quantum evolution in one map iteration is de- 
scribed by the unitary operator U: 

|$ = U\1>) = e ^ Tfi2 /2 e «.(§-*)'/2^ (3) 

where [9,h] = i, h = -id/89, and + 2tt)) = \ip(9)). 
Note that we have set K = 1. We study this map on the 
torus < 9 < 2ir, — tt < p < tt. The effective Planck 
constant is given by H e s = T. Indeed, if we consider the 
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FIG. 1: Poincare cross sections for the classical sawtooth map 
in the quasi-integrable regime at K — —0.5. We show 7 tra- 
jectories inside the integrable islands and a single trajectory 
filling the anomalously diffusive region. 

operator p = Th (p is the quantization of the classical 
rescaled action p) , we have 

[§,p]=T[9,n}=iT = ih eS . (4) 

The classical limit ft e ff — > is obtained by taking k — > 
oo and T — > 0, while keeping K — kT constant. We 
consider Hilbert spaces of dimension N — 2 nq , where n q 
is the number of qubits, and set T — 2ir/N. Therefore, 
h e g oc 1 /N — 1 /2 n « drops to zero exponentially with the 
number of qubits. 

The operator U can be written as the product of two 

operators, U k = e ik{ - § -^' 2 and U T = e- iTil ^ 2 . Since 
Uk is diagonal in the 9 representation, while Ut is diag- 
onal in the n representation, the most convenient way 
to simulate map © on a classical computer is based 
on the forward-backward fast Fourier transform between 
and n representations, and requires 0(N "log AT) op- 
erations per map iteration. The quantum computation 
takes advantage of the quantum Fourier transform and 
needs <3((logA) 2 ) one- and two-qubit gates to accom- 
plish the same task pH |2^ . More precisely, it needs 
2n q Hadamard gates and Sn 2 , — n q controlled-phase shift 
gates. Therefore, the resources required to the quantum 
computer to simulate the evolution of the sawtooth map 
are only logarithmic in the system size N, and there is 
an exponential speed up, as compared to the best known 
classical computation. 

Any experimental realization of a quantum computer 
has to face the problem of errors, which inevitably set 
limitations to the accuracy of the implemented algo- 
rithms. These errors can be due to unwanted couplings 
with the environment or to imperfections in the quan- 
tum hardware. In this paper, we limit ourselves to con- 
sider unitary errors, modeled by noisy gates. Such noise 



results from the imperfect control of the quantum com- 
puter. For instance, in a NMR quantum computer the 
logic gates on qubits are simulated by applying magnetic 
fields to the system. If the direction or the intensity of 
the fields are not correct, a slightly different gate is ap- 
plied, though it remains unitary. In ion-trap quantum 
processors, laser pulses are used to implement sequences 
of quantum gates |33| . Fluctuations in the duration of 
each pulse induce unitary errors, which accumulate dur- 
ing a quantum computation. 

As we have stated above, the implementation of 
the quantum algorithm for the sawtooth map requires 
controlled-phase shift and Hadamard gates |2lL |22| . We 
choose to perturb them as follows. Controlled-phase 
shift gates are diagonal in the computational basis and 
act non-trivially only on the four-dimensional Hilbert 
subspace spanned by two qubits. In this subspace, we 
write each controlled-phase shift gate as C = £C , where 
C is the ideal gate and the diagonal perturbation £ is 
given by £ = diag(e le °, e lei , e 162 , e" 3 ). Therefore, the 
unitary error operator £ introduces unwanted phases. 
The Hadamard gate can be seen as a rotation of the 
Bloch sphere through an angle 5 = tt about the axis 
uq = (sin 6*o cos 0o, sin sin cos 6>o), where 9q = tt/4 
and </>o = 0, so that u = (l/y/2, 0, 1/^/2). Since each 
one-qubit gate can be seen as a rotation about some 
axis u, unitary errors tilt the rotation angle: uq — > 
u = (sin 9 cos cj), sin 9 sin 0, cos 9), where 9 = 9 + Vx and 
4> = <p{) + i>2- We assume that the dephasing parame- 
ters €i,Vj (i = 1,...,4, j = 1,2) are randomly and uni- 
formly distributed in the interval [— e, e]. We also as- 
sume that the errors affecting different quantum gates 
are completely uncorrelated: every time we apply a noisy 
gate, the dephasing parameters randomly fluctuate in the 
(fixed) interval [— e, +e]. We note that the memoryless 
unitary error model has been wi dely investigated in the 
literature, see, e.g., Refs. pi I3I IM I37L l38| . 

We will compare the effect of noisy gates ("quantum 
errors" ) with that of randomly fluctuating perturbations 
in the system's parameters ("classical errors"). Here we 
choose to perturb the kicking strength k in Eq. J3J as fol- 
lows: at each map step, k is slightly changed by a small 
amount 6k(t), which is randomly chosen in the interval 
[-5k, dk]. Consequently, 5K(t) = TSk(t) E [SK,+SK], 
where SK = TSk. As we have discussed in the introduc- 
tion, this perturbation models, to some extent, the effect 
of round-off errors in classical computation. 

We will consider the following initial conditions: 

• A coherent Gaussian wave packet: 

N-l 2 

|^) G= ^^e^^ + <"--)^|n), (5) 

n=0 

where (9o,no) is the center of the wave packet 
((9) = 9q, (n) = uq), A a normalization con- 
stant, and er 2 = (An) 2 = ((n — (h)) 2 ) the variance 
in the momentum representation [39j. We choose 
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FIG. 2: Scattering circuit. The top line denotes a single ancil- 
lary qubit, the bottom line a set of n q qubits, H the Hadamard 
gate, and W a unitary transformation. 



a 2 = N / (2ttL) in order to obtain an equal value for 
the variances in p and in 9, namely A9 Ap = h e ff, 
with A9 = Ap = \/Htf£. The wave vector (J2J is 
the closest quantum analog of a classical probabil- 
ity density, localized in a small region of the phase 
space, centered in (9q,pq) and of width a. We point 
out that, as shown in Ref. |4(j, it is possible to 
prepare efficiently a coherent state on a quantum 
computer. 

• A random wave vector \^Po)r = Sn=i c n\n), where 
the coefficients c„ have amplitudes of the order of 
1 / s/N (to assure the normalization of the wave vec- 
tor) and random phases. This state has no classical 
analogue. 

The fidelity of quantum motion can be efficiently eval- 
uated on a quantum computer, as discussed in Ref. |26j |. 
Here we show an alternative method, based on the scat- 
tering circuit drawn in Fig. |21 0, |42] • This circuit has 
various important applications in quantum computation, 
including quantum state tomography and quantum spec- 
troscopy |42(. It ends up with a polarization measure- 
ment of just the ancillary qubit. We measure o~ z or a y 
and the average values of these observables are 

(a z ) = Re[Tr(Wp)\, (a v ) = Im[Tr(VFp)], (6) 

where (o~ z ) and (o~ y ) are the expectation values of the 
Pauli spin operators a z and a y for the ancillary qubit, 
and W is a unitary operator acting on n q qubits, ini- 
tially prepared in the state p (see Fig. |2J). These two 
expectation values can be obtained (up to statistical er- 
rors) if one runs several times the scattering circuit. If 
we set p — |V'o)('0o| and W = (Z7*)^ U*, it is easy to see 
that 

f(t) = |<Mtr*)t (jt Ml 2 = lTr{ w p )f = (a,) 2 + K) 2 . 

(7) 

For this reason, provided that the quantum algorithm im- 
plementing U is efficient, as it is the case for the quantum 
sawtooth map, the fidelity can be efficiently computed by 
means of the circuit shown in Fig. |2j 



III. RESULTS AND DISCUSSION 

Hereafter we will call f c (t) and f q (t) the fidelity decays 
induced by classical or quantum errors, respectively. 



Let us first consider the fidelity decay f c (t), obtained 
under fluctuating perturbations in the parameter k of 
the sawtooth map. We will show that, under this type of 
perturbation, the fidelity decay exhibits a marked depen- 
dence on the simulated dynamics. In particular, qualita- 
tively different behaviors are observed depending on the 
chaotic or non-chaotic motion. 

We first consider the quasi-integrable regime —4 < 
K < 0. In this case the sawtooth map behaves, inside 
the main integrable island with fixed point (9,p) = (n, 0) 
(see Fig. ^| as a harmonic oscillator, with characteris- 
tic frequency v K — ui K /2ir = \J — K/2n. Therefore, in 
the semiclassical regime the quantum motion of coher- 
ent wave packets residing inside integrable islands closely 
follows the harmonic evolution of the corresponding clas- 
sical trajectories. In the central island this motion has 
period T = 2n/^/—K, while in the outer islands the pe- 
riod is multiplied by a factor which depends on the order 
of the corresponding resonances (for example, the two 
upper islands in Fig. ^ correspond to a second order res- 
onance, and inside them the period is doubled). 

Since the chosen perturbation affects the parameter 
K, the fidelity f c (t) is obtained as the overlap of two 
wave packets which move inside an integrable island with 
slightly different frequencies. In this case, we know 
that for a static perturbation [8K (t) = 8K] the centers 
of the two wave packets separate ballistically (linearly 
in time) and a very fast decay of quantum fidelity is 
expected as far as the distance between the centers of 
the two packets becomes larger than their width a. The 
type of decay is related to the shape of the initial wave 
packet. In particular, for a Gaussian wave packet a Gaus- 



sian decay is expected. If 8v K = v K _ 



denotes 



the frequency separation between perturbed and unper- 
turbed motion, the Gaussian decay takes place after a 
time t s cx o~/5v K . 

In this paper, we consider the case of a randomly fluc- 
tuating perturbation 8K(t) € [— SK, SK]. Therefore, the 
frequency v K+6K{i) of a classical trajectory following the 
perturbed dynamics is not constant. The relative dis- 
placement of this orbit with respect to the one described 
by the unperturbed dynamics (with a frequency v K ) is 
approximately given by a Brownian motion. The sepa- 
ration between the two orbits is proportional to the fre- 
quency difference 8v K . In this case the fidelity decay is 
again Gaussian, but in general it shows large random 
fluctuations from the Gaussian profile (see for example 
the upper curve in Fig. |3J), which depend on the noise 
realization. Moreover, the distance between the centers 
of the two wave packets grows cx \/8v K t, and therefore 
the Gaussian decay starts after a time scale t s cx a 2 /8v K . 

Moreover, the fidelity decay depends not only on the 
shape of the initial state, but also on its position. Indeed, 
inside any integrable island the frequency's perturbation 
8v K = v K+SK — v K w A J^_ K is independent of the po- 
sition of the wave packet in phase space. Since larger 
orbits imply a larger velocity, and consequently a larger 
relative ballistic motion of the two wave packets, the fi- 
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FIG. 3: Fidelity decay for the quantum sawtooth map with 
n q — 12 qubits, in the presence of a classical fluctuating per- 
turbation in the k parameter. The initial condition is a Gaus- 
sian wave packet centered in (#o,Po) = (1,0). The upper 
curve shows the behavior in the quasi-integrable regime K = 
—0.5, with maximum perturbation strength 8K — 4 x 10~ 3 ; 
the lower one is obtained by simulating the map in the chaotic 
regime K — 0.5, with SK = 2 x 10~ 4 . In the inset we plot 
the same curves in a graph showing — log(/ c ) versus time. 
The straight lines correspond to exponential fidelity decay 
(— log/ c oc t, upper line) and Gaussian decay (— log/ c oc t 2 , 
lower line). 

delity drops faster when we move far from the center of 
the integrable islands. This is confirmed by our numeri- 
cal data (not shown here). 

In the chaotic regime, the fidelity f c (t) always decays 
exponentially, and an example of such decay is given in 
Fig- El For small perturbations, in the chaotic regime the 
decay rate T oc (5K) 2 , as predicted by the Fermi golden 
rule. However, if the perturbation is strong enough, the 
fidelity decay follows a semiclassical regime, in which the 
decay rate is perturbation independent and equal to the 
Lyapunov exponent of the underlying classical dynamics 
(see inset of Fig. 0J. The condition to observe the Lya- 
punov decay is that the perturbation is quantally strong, 
namely it couples many levels (8k > 1), but classically 
weak (5k <§C k). 

To summarize, the fidelity decay induced by classical 
perturbations strongly depends on the dynamical regime, 
chaotic or integrable. The two qualitatively different be- 
haviors (exponential or Gaussian decay) are shown in 
Fig. El Notice also that the regular dynamics turns out 
to be much more stable than the chaotic one (to rep- 
resent both cases on the same figure, the perturbation 
value chosen in the chaotic case is 20 times smaller than 
the one chosen in the integrable case). 

We now analyze the fidelity behavior in the presence of 
natural errors for quantum computation, namely random 
unitary perturbations of amplitude e on quantum gates, 
following the noise model described in Sec. [H] 



FIG. 4: Fidelity decay for noisy gates in the sawtooth map 
with K = 0.1, n q = 12. From right to left: e = 1.5 x 10~ 2 , 3 x 
10~ 2 , 4 x 10~ 2 , 5 x 10~ 2 , 6 x 10" 2 , 7.5 x 10" 2 , 10" 1 , 1.5 x 
10 _1 . Inset: fidelity decay for uncorrelated perturbations in 
the parameter k. From right to left 5K = TSk — 3 x 10~ 3 , 5 x 
10" 3 , 7.5 x 10~ 3 , 10" 2 , 1.5 x 10" 2 , 3 x 10" 2 , 5 x 10" 2 . In 
both graphs, data are averaged over 50 initial Gaussian wave 
packets. The two dashed lines show the Lyapunov exponential 
decay: f(t) = e~ At , where A ~ 0.315 is the classical Lyapunov 
exponent corresponding to K = 0.1. 
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FIG. 5: Characteristic time scale tf for the fidelity decay, 
determined by the condition f(tf) — 0.9, in the sawtooth 
map at K = 5, for the case of random noise errors in quan- 
tum gates. The data are obtained for different perturbation 
strengths e and number of qubits: n q — 4 (empty circles), 5 
(filled circles), 6 (empty squares), 7 (filled squares), 8 (empty 
triangles), 9 (filled triangles), and 10 (filled diamonds). The 
straight line shows the dependence tf ~ 0.126/e 2 ?! 2 ,, corre- 
sponding to the exponential fidelity decay JHJ, with C ps 0.28. 
The initial state is in all cases a Gaussian wave packet and 
data are averaged over 50 noise realization. 
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As sown in Figs. 1 1151 in the chaotic regime the fidelity 
f q (t) drops exponentially, with a rate T oc e 2 n 2 Q. This 
decay can be understood from the Fermi golden rule: 
each noisy gate transfers a probability of order e 2 from 
the ideal unperturbed state to other states. Due to the 
fact that perturbations acting on two different gates are 
completely uncorrelated, an exponential decay rate pro- 
portional to e 2 and to the number of gates n g = 3n 2 + n q 
required to implement one step of the sawtooth map is 
expected: 

f g (t) ~ e- Tt ~ e - c ^\ (8) 

where C ~ 0.28 is a constant which we have computed 
from our numerical data. We have determined the char- 
acteristic time scale tf for fidelity decay from the con- 
dition f q (tf) = A = 0.9 (note that the value chosen for 
A is not crucial). Our numerical calculations, shown in 
Fig. clearly demonstrate that: 

e n q 

in agreement with JS}. 

The fidelity decay in the chaotic regime always follows 
the exponential behavior predicted by the Fermi golden 
rule. Therefore, in contrast with the case of classical 
errors, there is no saturation of the decay rate to the 
largest Lyapunov exponent of the system (see Fig.^J. 

This result can be understood from the non-locality of 
quantum errors: each noisy gate can make direct transfer 
of probability on a large distance in phase space. This is a 
consequence of the binary encoding of the discretized an- 
gle and momentum variables. For instance, we represent 
the momentum eigenstates \n) (—N/2 < n < N/2) in the 
computational basis as \a Uq ■ ■ ■ 0.20-1) , where Oj € {0, 1} 
and n = -N/2 + N a j 2 ~ j ■ If we take ' sa ^ 
n q = 6 qubits (N = 2 6 = 64), the state 1 000000) cor- 
responds to \n = -32) (p = -tt), 1000001) to \n = -31) 
(p = -7T + 2tt(1/2 6 )), and so on until |111111), corre- 
sponding to 1 7i = 31) (p = — 7r + 27r(63/2 6 )). Let us con- 
sider the simplest quantum error, the bit flip: if we flip 
the less significant qubit (a\ = <-> 1), we exchange \n) 
with |n+l) (mod TV), while, if we flip the most significant 
qubit (a nq = <-> 1), we exchange |n) with |n + 32) (mod 
N). It is clear that this latter error transfers probability 
very far in phase space and cannot be reproduced by clas- 
sical local errors. Therefore, no semiclassical regime for 
the fidelity decay is expected with quantum errors. In 
particular, the non locality of perturbations makes the 
fidelity insensitive to the rate of local exponential insta- 
bility, given by the Lyapunov exponent. 

The most striking feature of the fidelity decay induced 
by quantum errors is that it is substantially independent 
of the chaotic or non-chaotic nature of the underlying 
classical dynamics. An example of this behavior is shown 
in Fig. and strongly contrasts with what obtained by 
perturbing the system's parameters (see Fig. [3}. In par- 
ticular, the fidelity decay for integrable dynamics is ex- 
ponential, as shown in Fig.|SJ If we start from a Gaussian 




FIG. 6: Fidelity decay for the quantum sawtooth map sim- 
ulated with n q = 12 qubits, in the presence of uncorrelated 
unitary quantum errors with maximum perturbation strength 
e — 10~ 2 . As initial condition we consider a Gaussian wave 
packet peaked in (#o,f»o) = (1,0). The upper curve shows the 
behavior in the quasi-integrable regime K = —0.5, while the 
lower one is obtained by simulating the map in the chaotic re- 
gion K — 0.5. In the inset we plot the same curves, showing 
— log(/ 9 ) versus time. The solid line corresponds to exponen- 
tial fidelity decay, that is — log f q oc t. 

wave packet, integrable dynamics turns out to be a little 
more stable than chaotic dynamics: we numerically ob- 
tained a ratio of the decay rates in the chaotic and in the 
integrable case which oscillates between 1.15 and 1.4, for 
different values of n q between 5 and 16, and for various 
e ranging from 10~ 5 to 10 _1 . 

We stress that the smaller decay rate obtained when 
we evolve a Gaussian wave packet inside an integrable 
island is not due to the lack of exponential instability 
but simply to the fact that the dynamics preserves the 
coherence of the wave packet. This can be clearly seen 
from the data of Fig. [7| 

• In the chaotic regime K > (Lyapunov exponent 
A > 0), the fidelity decay rate is independent of the 
initial state (Gaussian packet or random state) and 
of the rate of exponential instability. Indeed, the 
decay rate is independent of K, while A depends on 
K. 

• In the quasi-integrable regime —4 < K < (Lya- 
punov exponent A = 0), only in the case in which 
we choose as initial state a Gaussian packet placed 
inside an integrable island we obtain a fidelity de- 
cay rate smaller than in the chaotic case. On the 
other hand, if we start from a random state of if we 
place the Gaussian wave packet inside the anoma- 
lously diffusive region, we obtain the same decay 
rate as in the chaotic case. 

From these results, we conclude that the decay rate does 
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FIG. 7: Dependence of the fidelity decay rate, induced by 
quantum uncorrelated unitary perturbations, on K, for n q — 
9, e = 10~ 2 . The dashed line separates the quasi-integrable 
region — 4 < < from the chaotic region K > 0. As 
initial condition we choose: (i) a Gaussian wave packet cen- 
tered in (8o,po) = (1,0) (circles) (note that for —4 < K < 
this packet is inside the main integrable island); (ii) a Gaus- 
sian packet centered in (#o,Po) = (0,0) (squares), that is re- 
siding in the diffusive region; (iii) a random wave function 
(diamonds). All data are obtained after averaging over 25 
different noise realizations. 



motion. The main result is that, while the fidelity decay 
under classical errors strongly depends on the dynamical 
nature of the system under investigation and on initial 
conditions, quantum errors act in a way essentially inde- 
pendent of the system's dynamics. This practical insensi- 
tivity to the dynamics is eventually a consequence of the 
intrinsic non locality of the errors that naturally affect 
the quantum computation. As a consequence, the rich 
variety of behaviors found in the study of the stability 
of quantum motion under perturbations of the s yste m's 
Hamiltonianjj], B II S S II 1 S E3 El El M III 
ITU liH. IrH Il8l| has no correspondence in the fidelity of 
quantum computation under its natural perturbations. 
The stability of quantum computation is essentially in- 
dependent of the chaotic or integrable behavior of the 
simulated dynamics. This conclusion is simply based on 
the non locality of quantum errors and therefore we ex- 
pect that it remains valid also in the case of non-unitary 
quantum noise and/or when errors, correlated or mem- 
oryless, act not only on the qubits on which we apply a 
quantum gate but on all the qubits that constitute the 
quantum computer. 



not depend on the value of the Lyapunov exponent. In 
short, the decay of the fidelity due to noisy gates is inde- 
pendent of the presence or lack of exponential instability 
|45| . We point out that we have checked that this state- 
ment remains valid also for static errors, like in the case 
in which the dephasing parameters ti , Vj appearing in our 
noise model are time- independent. 

IV. CONCLUSIONS 

In this paper, we have compared the effects of clas- 
sical and quantum errors on the stability of quantum 
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